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Abstract 

In the framework of HQET inclusive decays of b-flavored hadrons can be 
handled as decays of the free b-quark approximately. We analyse the correc- 
tion to this approximation for the inclusive decays into a polarized J/ip, the 
correction is characterized by two matrix elements defined in HQET. For the 
J/ip we use NRQCD to parameterize the formation of the cc pair into J/iJj. 
Numerically the correction is remarkably large, it can be at the level of 30% 
and it is mainly due to the Fermi-motion of the 6-quark. With this correc- 
tion we give a new determination for combinations of two NRQCD matrix 
elements. 
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Inclusive decay of a b-flavored hadron into a charmonium is an interesting process in the 
aspect that two effective theories derived from QCD are used to factorize nonperturbative 
effects. For the charmonium, because the relative velocity v c between the c- and c-quark 
inside the charmonium is small in the rest-frame, one can use NRQCD |T[to handle the 
nonperturbative effect in the formation of the free cc pair into the charmonium, while for the 
initial hadron one can employ an expansion in the inverse of m&-the b-quark mass, because 
rrib is very large in comparison with the nonperturbative energy-scale A. The expansion 
can be done with the heavy quark effective theory(HQET) 0(for a review, see e.g. ||). 
At the leading order the decay can be considered as the decay of the free b-quark, hence 
the decay width is the same for all initial hadrons. Corrections from higher orders can be 
systematically added, they will depend on the type of initial hadrons. It is the purpose of 
this work to study the corrections from the next-to-leading order, i.e., the effect of initial 
hadrons. 

It should be kept in mind that the factorization mentioned before is not well estab- 
lished in comparison with the factorization for totally inclusive decays and for inclusive 
semilepton-decay, in these decays one can use the operator product expansion(OPE) to fac- 
torize nonperturbative effects. With OPE the effect of the initial hadron in inclusive decays 
are studied, for example, in @||. For the process considered here OPE can not be used. 
However we can use the diagram expansion proposed in |J to perform the factorization and 
assume that the factorization still holds, i.e., the perturbatively calculated coefficients are 
free from infrared divergences. We will use the expansion and work at the leading order of 

We consider the process: 

H b (P)^J/^j(k,X) + X (1) 

where we denote the b-flavored hadron as Hf,, which contains a b-quark. The momenta are 
given in the brackets. A denotes the polarization of the J/ip, A = L is for the longitudinal 
polarization and A = T is for the transversal polarization. To analyze the decay three 
expansions are used, they are expansion in a s , in v c and in the inverse of m^. In |7],[| the 
decay is studied at the leading order of all expansions. At the leading order of v c and of 
the inverse of rrib the one-loop correction is calculated in || for the unpolarized J/ip. The 
corrections from the next-to-leading order of v c is analyzed in fill at the leading order of a s 
for the polarized J/ip, they can be very large. Because of lacking the detailed information 
about the formation of the cc pair into a J/if>, i.e., the precise values of NRQCD matrix 
elements, a detailed prediction seems impossible without help of some models. In this work 
we analyze the corrections from the next-to-leading order of the inverse of rrib to complete the 
analysis of the correction from the next-to-leading order of all expansions. By completing 
this analysis the last unknown effect in the theory for the process can be estimated, in which 
nonperturbative effect is represented by matrix elements defined in HQET and in NRQCD. 
In [11 1 models to account the effect are used to explain the J/^-spectrum. 
The effective weak Hamiltonian for the decay is: 

H * = % E {V cb V: q [\c {1] ^)cY^ - 7fe)cyy„(l - 7 5 )& 

+q 8] (/i)cTV(l - 7 5 )cgT a 7 ,(l - 75)6]}- (2) 
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We neglected the contributions of QCD penguin operators in if e fr- T a (a = 1, • • • 8) is SU(3) 
color-matrix. The coefficients C[i] and C[%] are related to the usual C± by 



C w (/i) = 2C+( M ) - C_(Ai), C7 M (Ai) = C + (m) + C-0*)- (3) 

With the effective Hamiltonia a cc pair can not only be produced in the color-singlet state, 
but also in the color-octet state. Although a cc pair in the color-octet state can be trans- 
mitted into a J/ij) at higher orders in v c than a cc pair in the color-singlet state does, its 
contributions to the decay rate are important, because the cc pair is more likely in the color- 
octet state than in the color-singlet state by the fact C[i](ra&) : C[s](m b ) ~ 0.18 : 1. We will 
take the color-octet states into account and write the decay width as 

rA(J/^) = ri 1) + ri 8) , (4) 

where the index 1 or 8 stands for color singlet contributions or for color-octet contributions, 
respectively. With the H e ff and in the approach of the diagram expansion at the leading 
order in a s , the contributions to the decay width can be written: 



X 



■ (015(0)7^(1 - T5 )c(0)| J/iP + X)(J/tP + X\c(x)Y(l - 7s)c(a;)|0) 

• (# 6 |6(0) 7m (1 - 75)7 • qiAl - K)b(x)\H b ), (5) 

where we used \V c d\ 2 + \ V CS \ 2 ~ 1. The integral J dT is 

In the above equations we neglect the mass of light quarks and take the nonrelativistic 
normalization for Hf,- and J/^-state. Similarly the color-octet contribution is 



X 

b„Mt 



■ (0|5(0)TV(1 - T5)c(0)| J/ip + X){J/^j + X\c(x)T b Y(l ~ 7s)c(x)|0) 
• (F 6 |6(0)T a 7 M (l - 75)7 • qT b l„(l - ls)b{x)\H b ). (7) 

In both contributions the dependence on the initial hadron appears in a matrix element, 
which is given by: 

Tijix) = {H b \bMh{x)\H h ) } (8) 

where the indices i,j stand for color- and Dirac indices and the average of the spin is implied 
if H b has a non-zero spin. This matrix element represents the nonperturbative effects related 
to the initial hadron, it can be expanded in the inverse of m b in the framework of HQET. The 
expansion is performed to expand the Dirac fields 5(0) and b(x) with the fields in HQET. We 
use the expansion proposed in |0| , which is slightly different than the standard approach in 
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0. However, the final results are same by using these two expansions, if one correctly takes 
the normalization of the state in HQET into account. The Dirac field can be expanded as: 

b(x) = {x + J_ i7 . Dt + _L t , . Dl . Dt _ _i f(7 . DT y} h(x) + 0( i f) 

+ (terms for anti — quark), (9) 

where is the velocity of H b , is the covariant derivative + igG ti (x). Dj. is defined as 

D»= D» -v^v D. (10) 
h(x) is the field of HQET and its equation of motion reads: 

{ iB . c __i. (7 . CT)2 } ftW = 0+O (-l). (id 

With the expansion in Eq.(9) the operator N for the number of b-quarks is 

N = f d 3 xh(x)h(x). (12) 



Using the expansion in Eq.(9) the matrix element can be written as 



1 



Tij{x) = e- m » v - x {{H b \hMhi{x)\H b ) + — (# 6 |^(0)(i 7 ■ D T h(x)) t - t 1 ■ Drh^h^x) \H b ) 



+ ^(H b \hj(0)(v ■ D 1 ■ D T h(x)) t + v-D 1 - D T hj(ti)hi{x) + 7 ■ Brfy(0)(7 • D T h(x))i 

-i^(0)(( 7 • £>t) 2 /i(x)), - \'^D^h{Q)h i {x)\H b )} + O(ij). (13) 

The dominant x-dependence of Tij(x) is exp(— im b v ■ x), while the ^-dependence of matrix 
elements with h(x) is controlled by the energy scale A, which is much smaller than m b , we 
can expand this x-dependence in 1. In this expansion we obtain a tower of local operators 
O n . The equation of motion Eq.(ll) should be used to eliminate redundant operators. The 
contribution from a given operator O n to the decay width will be suppressed by the power 
A^ dn ~ 3 \ where d n is the dimeneion of the operator O n . For example, the first matrix element 
in Eq.(13) is expanded as 

(tf 6 |£ i (0)^(aO| J H 6 ) = (^1^(0)^(0)1^) + x^{H b \h s (0)& t h i (0)\H b ) 

+^x v (H b \h j (0)d^d v h i (0)\H b ) + ■■■. (14) 

With the Lorentz covariance and properties of the field h(x), the matrix element can be 
written 

(Hblhj^hiix^Ht,) = ^(1 + 7- v) l3 {H b \hh\H h ) + ^vx(l + 1 - v^H^hv ■ dh\H b ) 

+^(1 + 7 ■ v) ijXtl x v [g^{H b \hd^h\H b ) 

-v^v v (H b \h(d^ - 3{v ■ dfh\H b )} + ■■■. (15) 
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Because of Eq.(12) the first matrix element in the r.h.s. of the above equation is 1: 

(H b \hh\H b ) = 1. (16) 

This results in that the decay of H b at the leading order of m^ 1 is the decay of the free 
b-quark. Similarly we can expand the other matrix elements in Eq.(13) and we obtain: 
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Ty{x) = e-™*~{-{! + 7 " «)« + — [-(1 + 7 



12 v ' , J 2m fe L 12 
+ T^(7m - V? ■ v)]{H b \h(l ■ D T fh\H b ) + ^x v {g, 

■(1 + 7 • v) i:j (H b \hD 2 T h\H b ) - ^{H b \hin ■ D T fh\H b )} + ■■■ (17) 

where we have replaced the derivative <9 M with the covariant one and the equation of 
motion was used to eliminate the redundant operators. The ■ • • stand for contributions from 
operators with the dimension higher than 5, which will lead to contributions at the order 
higher than A 2 . 

For the chamonium J /if} we will work at the leading order in v c and we have Mj/^ = 2m c . 
To expand the Dirac field c(x) with the field of NRQCD, we make a Lorentz boost to 
transform the relevant matrix element into the J / if>-iest frame and then make the expansion. 
We obtain for the color-singlet part: 

£(0|c(0)7"(l - 7s)c(0)| J/^ + X)(J/if> + X\c(x)Y(l ~ 7s)c(x)|0) 
x 

= ^^e^XWiXW^Oi'^'S,)) + 0(v 2 ), (18) 

where the factor Mj/^/k is due to the boost. The definition of the matrix element 
(Ol^fSi)) can be found in §. 

With Eq.(17) and Eq.(18) we can perform the x-intergration in Eq.(5). After the inte- 
gration we obtain a ^-function for the momentum conservation, and also derivatives of the 
5-function like 

JL(2Trf8\m b v-q-k). (19) 

Such derivatives can be eliminated by partial integrations when the integration of the phase- 
space is performed. For example, the contribution related the above derivative can be 
written: 

d 4 qS(q 2 )f(q)^5\m b v - q - k) = -(^(q 2 )f(q) + 5(q 2 )^f(q))\ q=mb ^ k , (20) 
oq^, oq^ 

where 8'{x) = 4-8(x). Performing these integrations the final results for the color-singlet 
contribution can be obtained. To present our results we write: 

(H b \hD 2 h\H b ) = 

(if 6 |% • D T ) 2 h\H b ) = (H b \hD 2 T h\H b ) + (H b \h~gG^h\H b ) + ■■■ 

.2 , ,.2 



/4+K + ---, (21) 



where • • • stand for higher-order contributions. The color-singlet contribution can be 
written as 

r?) = >^^ {0 , rs ^ . { Fx(y) + a hm + | GaM } (22) 

where y = 4ml /ml an< l the functions are: 

F L (y) = (l-y) 2 , 
F r ( 2/ )=y(l- 2/ ) 2 , 

G L (t/) = i(7-% + 7y 2 -10t/ 3 ), 
o 

G T (y) = ^(ll-26y + 15y 2 ), 
^L(y) = ^(-53 + %-y 2 + 26y 3 ), 

H T (y) = \y (-11 + l%-27y 2 ). (23) 

From the above results the effect of the initial hadron is represented by two matrix elements 
defined in HQET. It should be noted that the effect is not suppressed by the power of m^ 1 , 
but rather by the power of (m 2 — Ami)' 1 . This can be seen by comparing the coefficients 
given above. 

For the color-octet contribution we expand the the color-octet matrix element related to 
J ftp with NRQCD fields. The result at the leading order of v c is: 

5>|c(0)TV(l - 75)c(0)]J/^ + X){J/*I> + X\c(x)T b Y(l ~ ^)c(x)\0) 
x 

~ Jk - m MTT-(O J 8 '*eS )) +^(A)(^(A))*(0 8 J/ *( 3 5 1 )> 



fco 24 M 



k »k u „ 1 



H-9T + ^(A)(^(A))*)- i (0 8 ^( 3 P 1 ))} + 0{vt). (24) 



At this order, the color-octet cc pair with the quantum number 1 Sq, 3 Si and 3 Pi can form 
a J/ijj through emission or absorption of soft gluons, the probability is at order of v\. 
Perform similar calculations as for the color-singlet contribution we obtain the color-octet 
contribution: 

where the functions are: 

*f(v) = (l - y) 2 «O 8 J/ Ys )> + (O^CS,))) + 2y(l - yf^iOPfPi)), 

m c 

F?\y) = (i - vmoi'+fSo)) + i/<o 8 J/ *( s Si)» + (i - y)\i + y^io^epj), 

m c 
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4 8) (?/) = g((7 - 12y + I5y 2 - 10y 3 )(O J s ^CS )) + (7 - Ay + 7y 2 - 10y 3 )(O J s / ^ 3 S 1 )) 
+2(-i - Uy + lBy 2 )-^^^))), 

G?(y) = \({7- 12y + I5y 2 - 10y 3 ) (O^ C So)) + y{U - 26y + 15y 2 ) (0 J 8 /lp ( 3 S 1 )) 
6 

+ (7 - 17y + 5y 2 + 5y 3 )^(0 8 W ( 3 P 1 ))), 

H?\y) = + 12y - 9y 2 + 26y 3 )(O 8 J/ ^( 1 5 )) + (53 + Ay - y 2 + 26 2 / 3 )<0 8 7 ^( 3 S' 1 )> 

o 

+2y(l + 2y-27y 2 )^(0^( 3 P 1 ))), 
m A c 

H [ t\v) = ~((-53 + 12y - 9y 2 + 26y 3 ) (O^ C S )) + y{-U + Uy - 27y 2 )(0^ \ 3 S l )) 
6 

+ ( _ 53 + l7y -llj/ 2 - z/ 3)i_ ( o 8 / ^( 3 P 1 ))). (26) 

m c 

With the results given above we complete our analysis for the effect of the initial hadron. 
As mentioned before, because the NRQCD matrix elements, especially the color-octet matrix 
elements, are not known precisely and also there are possibly large relativistic corrections 
for charmonium, we can not give a precise prediction for the decay width. However, the 
effect of the initial hadron can be estimated numerically. For this we take m& = 4.8GeV and 
m c = 1.5GeV. We obtain for the color singlet contribution: 

rf = m > G ^ V *\ (0^(3^)) . jo.15 - 0.63 
4327rm c [ 

if = ^^^'V / 7 ^ 3 ^)) ■ {o-37 - 8.34^ + .9 8 £| \ . (27) 

It is remarkable that the coefficient in the front of /j, 2 is very large for 1^ in comparison 
with others. This indicates that the effect of the initial hadron may be substantial. Taking 
Hb = B as an example, the parameter ^ 2 g {B) can be determined by the mass splitting 
between B and B* |Tj|, it gives ^ 2 (B) fa 0.36GeV 2 , while fJ> 2 {B) can be determined from 
the QCD sum rules WM and from an analysis of spectroscopy of heavy hadrons [I5|. These 
determinations give a numerical ranee from 0.3GeV 2 to 0.54GeV 2 for nl{B). There is also 
a constraint for fil(B), f^l(B) < ^ 2 (B) [fTBl . We take the value [i 2 {B) « fi 2 (B) to estimate 





+ 0.20-^4 


m\ 




A 


+ 0.98^-| 


m\ 


mi 



the effect of the initial B. With these values, the effect in 1^ is at 4% level, which is 



negligible, while the effect in 1^ is to reduce the width at 30% level in comparison with 
at the leading order, which is significant. Taking the same quark masses we obtain for 
the color-octet contribution: 

rf = ^Cwl^l 2 . {(0 .37 - 8.024 + o.674)<of^o)> 
2887rm c mi mi 

2 „2 

+ (0.15 - 0.63^ + 0.20H) (O^CSt)) 
mi mr 
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2 2 i 

+ (0.52 - 8.02^ + 0.23% — (0^( 3 Pi))} 

Traj; 

r ( 8) = rjglgl^ _ ^ /i )(0 ^ (ls o)) 

L 2887Tm c u m 2 m 2 ^ 8 v u;/ 

2 u 2 

+ (0.37 - 8.34^ + 0.98% (O^fSr)) 
ml mi 

2 2 i 

+ (0.29 - 0.30^ + 0.54^|) — (O^fPO)}. (28) 

mj. Tra- 
in the color-octet contribution there are also large numbers in the front of /j 2 in different 
production channels. This results in that the corrections are significant. We also take 
P-meson as an example and we find: 

r (8) (jB) w . {(0 .37 - 0.11)(0 8 ^(^o)) 

2887rm c 

+(0.15 - 0.007) (Oi'^SJ) + (0.52 - 0.12)^ ¥ (O 8 J/,/ '( 3 P 1 ))} 

m c 

r ,„ w . {(0 .37 - 0.11)<Q**(%)> 

2887rm c 

+ (0.37 - 0.11)(O^( 3 S' 1 )> + (0.29 - 0.004)^-(O 8 J/,/, ( 3 P 1 ))}, (29) 

m A c 

where the second number in (■ ■ •) is for the corrections. From these numbers one can see 
that the corrections are large, especially in the 1 Sq channel, which can be at 30% level. The 
large corrections arise in both color-singlet and color-octet contributions from terms with 
/x 2 , that indicates that the Fermi-motion of the 6-quark inside B affects the decay width 
substantially. 

If we collect all known results to compare the experimental result [TjJ 

Br[B -> J/tp + X) cxp = 0.80 ± 0.08%, 
Pr(P -»• ip' + X) exp = 0.34 ± 0.05%, (30) 

we have serious problems with the one-loop QCD correction and with the relativistic correc- 
tion. With the one- loop correction the color-singlet contribution becomes negative M. The 
relativistic correction is analyzed in [K| and it can have large effect, but the detailed size is 
unknown because there are four unknown matrix elements of NRQCD. If we neglect it we 
have: 



Br(B -f if, + X) = 0.073 • 10" 2 (Of ( 3 5i)) 



0.34 



-0.19(Of ( 3 5x)> + 0.33(OtCS )) + -^(Of ( 3 P)) 

m z c 

-{0.068 • 10~ 2 (Of ( 3 5i)) + 0.071 (OtCS )) 



+0.027(O?( 3 5 1 )) + ^(Of ( 3 P))}, (31) 
mi 
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where we use if) to denote J/ip or if)'. In the first two lines there are the results from the 



leading order in the -expansion. The numbers are from || and they contain one-loop 
QCD corrections. They are slightly modified due to the m^ 2 corrections in the semiletonic 
decay width. It should be noted that the color-singlet contribution is improved by adding 
certain contributions from higher orders and with the improvement it becomes positive. The 
results of this work are given in the last two lines. With the corrections the contributions 
from all channels become smaller than those without the corrections. If we take a larger 
value of /j%(B), the color-singlet contribution becomes negative. Among the four matrix 



elements in Eq.(31) the color-single one is calculated with potential models |18fl and with 
lattice QCD |l!|, whose value is fii 



(O^CSt)) = 1.16GeV 3 , 
(of( 3 Si)) = 0.76GeV 3 . (32) 

The matrix element (Of^So)) is determined by experiments at Tevatron and at Hera 
fl27H2~3"l , but the uncertainty can be large. We take the value: 

(O^CS^) = 1.06 • l(T 2 GeV 3 , 
(Of ( 3 Si)) = 0.44 • 10~ 2 GeV 3 . (33) 

The other two are not well determined, only certain combination of them is known with 
large uncertainty. With the experimental results in Eq.(30) we can determine a combination 
of them: 

(O^eSo)} + ^(O^fPO) = 2.4 ■ 10- 2 GeV 3 , 

(Of? So)) + ^-(Of CP,)) = 1.0 • 10~ 2 GeV 3 . (34) 



c 



Comparing the combinations determined without the corrections the change is significant, 
the combinations becomes 80% larger than those without the corrections and they are closer 



to those determined from other experiments [20-23|. However, the determination of the 
combinations should be regarded as a rough one because possibly large corrections from 
higher orders in v c are neglected. 

To summarize: We have analyzed the effect of the initial hadron in the process — > 
ip + X. With our results the three expansions used for the process are all completed at the 
next-to-leading order. The effect of the initial hadron can be at 30% level for B —>■ ip + X 
and it reduces the decay width. The effect is mainly due to the Fermi-motion of the 6-quark 
inside the initial hadron. Including this effect we have determined combinations of two 
NRQCD matrix elements, which are 80% larger than those without the effect and they are 
closer to those determined in other experiment. However, the results for the combinations 
should be taken with caution because neglected effects can be large. 
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